Photoinduced magnetism in the ferromagnetic semiconductors 
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We study the enhancement of the magnetic transition temperature T c due to incident light in fer- 
romagnetic semiconductors such as EuS. The photoexcited carriers mediate an extra ferromagnetic 
interaction due to the coupling with the localized magnetic moments. The Hamiltonian consists of 
a Heisenberg model for the localized moments and an interaction term between the photoexcited 
carriers and the localized moments. The model predicts a small enhancement of the transition 
temperature in semi-quantitative agreement with the experiments. 
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There is a class of materials which are ferromagnetic 
semiconductor s! 1 ' 2 ' 3 ' ' 5 such as the Europium chalco- 
genides EuX (X=S, O, Se, Te) and chalcogenide spinels 
MCr 2 Y 4 (M=Hg, Cd; Y=S,Se). The ferromagnetism 
in these materials are generally well described by the 
Heisenberg model. In this paper, we are concerned with 
the effect of incident light on the magnetism, in par- 
ticular, on the ferromagnetic transition temperature T c . 
The basic mechanism of the photoenhanced magnetism 
is that light excites electron-hole pairs, which mediate 
an extra magnetic interaction between the local mag- 
netic moments, on top of any pre-existing magnetic in- 
teraction, thereby enhancing the T c . The related phe- 
nomena, where ferromagnetism occurs only if light is 
present but otherwise the system is paramagnetic, has 
been studied in connection with the dilute magnetic 
semiconductors^^^i 3 .^^ 

Here we study the photoinduced magnetism in the fer- 
romagnetic semiconductors by solving a model Hamil- 
tonian within a mean-field approach by constructing a 
BCS-type wave function for the electron-hole pairs^ 
similar to our analysis developed earlier for the di- 
lute magnetic semiconductors^ Our work provides a 
semi-quantitative explanation for the enhancement of T c , 
which has been observed to be as large as 0A6K for 
EuS£»Z»& under intense laser radiation. 

The ferromagnetic semiconductor is modeled by 
isotropic conduction and valence bands interacting with 
fixed Eu magnetic moments originating from the f elec- 
trons. The bands are spin-split through an exchange 
interaction with the magnetic ions as indicated in Fig. 
[TJ The Hamiltonian consists of a Heisenberg interaction 
term between the fixed Eu magnetic moments and a term 
describing the itinerant band electrons and their interac- 
tion with light. We thus write 




H = — ^ J M Si-Sj + Ti, 



ell i 



(1) 



where the first term is the Heisenberg term and < ij > 
denotes summation over distinct pairs of nearest neigh- 
bours. Each Eu atom has twelve nearest neighbor in the 
fee lattice of EuS. The second term H^ h describes elec- 
trons and holes making up the conduction and valence 



FIG. 1: A schematic picture of the photoinduced ferromag- 
netism in EuS. Incident light creates electron-hole pairs that 
mediate an extra ferromagnetic exchange interaction between 
the Eu moments. The net result is an enhancement of the 
critical temperature T c by a fraction of a degree K for typical 
parameters. 



bands and their interaction with light 

Hjf, = Kkin +H C + (2) 
The kinetic energy term is given by 

H kin = J2( E kA a c ka + E h k J ka d ka ), (3) 

ha 

where c\ a , Cka and d\ a , dka are the field operators for the 
electrons and holes, respectively, with k being the Bloch 
momentum and a, the spin index. Throughout this pa- 
per, electrons refer to the photoexcited electrons in the 
conduction bands and holes refer to the lack of an elec- 
tron in the valence bands. Since the photon has zero spin 
and negligible momentum, absorption of light involves a 
vertical transition of an electron in the valence bands into 
the conduction bands with the spin of the electron con- 
served thus resulting in the creation of an electron-hole 
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pair described by cl cr d_ k _ a . Without incident light, we 
have the "vacuum state" with no electrons nor holes. 

The single-particle energies in Eq. [3] include the Zee- 
man splitting due to the electron interaction with the 
localized Eu spins, so that 

E ka = E 9+ ± Je <S> (4) 

2m e 

and 

n 2 k 2 

El = — ±J h <S>, (5) 

where ± describes the Zeeman splitting of the up-spin 
and the down-spin states. Here J e and Jh are the elec- 
tron and hole exchange interactions with the localized Eu 
moments, < M >— gps < S > is the average moment of 
a Eu atom, m e and are the electron and hole masses, 
respectively, and E g is the band gap. 

The Coulomb interaction part consists of three terms 

H c = H ee + T~thh + T~teh, (6) 

where the electron-electron interaction term is 

Hee = 2 ^ ^ C t+q,a C l'-q,a' C k'^C k a, ( ? ) 
kk' qaa' 

with a similar term for the hole-hole interaction, while 
the electron-hole interaction is given by the expression 

U eh = - ^ V g cl +q a dl^ q a ,d k , a ,c krT . (8) 

kk' qaa' 

Finally, the last term in Eq. (|2|) which describes 
the coupling to the external electromagnetic field of fre- 
quency oj, in the "rotating wave approximation" is given 
as 

n L (t) = A^p fe [(4 T dL fci +4A T ) e " iwt + ( 9 ) 

k 

As mentioned earlier, the light induces a vertical tran- 
sition with no net spin change and Xp k is the cou- 
pling strength for the creation of such an electron-hole 
pair. The coupling strength is given as \p k = eE < 
ipvk\e ■ pl^Pck > /muj where ipek, ipvk are conduction and 
the valence state wave functions, E is the strength of the 
electric field and e is the electric polarization vector and p 
is the momentum operator. The matrix element may be 
computed once the wave functions are known, but here, 
for simplicity, we neglect the momentum dependence of 
the matrix elements, so that p k = 1 and the parameter A 
indicates the strength of the interaction. For strong laser 
radiation available in the laboratory, a broadly estimated 
value used by earlier author o 16 i 18 is A « 0.1 eV, which 
we use here in our numerical calculations. 

The time-dependence of Eq. §§§ can be eliminate d 14 ^ 8 
by an appropriate unitary transformation e~ lS TLe lS 



where S = -^J2ka( c la c ka + d\ a d kc! ) and the Hamil- 
tonian becomes 

a = — ^ jhiSi.Sj 

<i,j> 

+ ^(Ka - f^/Ma^a + {E& - hw/2)d\ a d ka 

+H C + A^{p fe (4/_ fci + 4/_ feT ) + h.c.}. (10) 

Thus the time has been eliminated and the transformed 
Hamiltonian may be interpreted as a "quasi free en- 
ergy" with the chemical potential for the creation of an 
electron- hole pair given by fi — hw. One can thus write 
the transformed Hamiltonian in the form Ti — H.QE—pN, 

where N = 2\lka c \a Ck ° 1S ^ ne number of electron-hole 
pairs. 

The Coulomb interaction term appearing in the Hamil- 
tonian ([7]) and ([5]) is given by V q = Aire 2 / (eq 2 ), but in our 
discussions below, we replace it, for the sake of simplic- 
ity, by a contact interaction term V(r) — VoS(r) without 
changing the qualitative physics. Its typical magnitude 
in semiconductors is Vq ~ 20 eV.A 3 ^ Typical value 
for the exchange parameter is: Jm = 2.19 x 10~ 5 eV, 
which is calculated using the expression for the criti- 
cal temperature for the Heisenberg Hamiltonian /cbT c = 
(v/3)J M S(S + l). For EuS, T c = 16K, S = 7/2, and the 
number of nearest neighbours v — 12. The magnitude 
of J = J e — Jh = —0.036 eV as obtained from estimat- 
ing J e and Jh from the splitting of the spin-up and the 
spin-down bands in the band structure of EuS. 1 

We compute the free energy of the system by using a 
variational wave function similar to the Bardeen-Cooper- 
Schricffcr (BCS) wave function with zero net momentum 
and spin. The variational wave function consists of a 
superposition of the electron-hole pair a 19 i 20 generated by 
the coupling to light 

I* >=Y[(u k a+VkacU_ k _ a )\0>, (11) 

ka 

where |0 > is the intrinsic vacuum state (filled valence 
and empty conduction bands), v k<T is the probability am- 
plitude for the creation of an electron-hole pair with the 
electron momentum and spin of ka and a hole with the 
opposite momentum and spin, while u k<7 is the proba- 
bility amplitude for the absence of such pair excitation. 
The normalization condition is, as usual, u ka . + v ka = 1, 
and the coefficients u ka and v ka may be chosen as real. 
In the present case the localized moments break the sym- 
metry between the spin up and down electrons leading 
to a net spin polarization of the electron and holes. This 
in turn affects the net magnetic moment, which then af- 
fects the energies of the particles through the Zeeman 
term, so that all quantities have to be determined in a 
self-consistent manner. 

We now calculate the expectation value of the "free 
energy" from Eq. (|10[) , which may be rewritten in terms 



3 



of the gap parameter to yields 

F = + il a /2)v% a - (A a u kr7 v ka /2 + \u ka v k<J )] 

(12) 



-3 M <S> 2 xAf, 



where 



and 



A CT = 2Vq ^2 u k'aV k 'a + 2A 



(13) 



(14) 



Here, v is the number of nearest neighbour (12 for the fee 
lattice) and the factor of 1/2 in the last term of Eq. [12] 
appears to avoid double counting. TV is the total number 
of sites in the crystal and £, ka is the energy of the electron- 
hole pair excitation with a spin-up electron and a spin- 
down hole measured with respect to the photon energy 



E 



= E n 



El 

h 2 k 2 

2m* 



±J<S> -hid, 



(15) 



where J — J e — Jh and the effective mass of the electron- 



hole pair m* 



Minimizing the free energy 



with respect to v ka , dFo/dv ka = 0, and using the normal- 



ization conditions 



= 1, we get the occupation 



amplitudes u ka and v ka , which read 



4 = ~[1 



e 



and 



2 [ y/Ai + (t ka + ny 



(16) 



(17) 



We solve the equations (fl"6|) and (| 1 T[) numerically which 
yield the results for u ka and v ka , from which the free 
energy Fq is computed using Eq. (|12|) . The typical pair- 
occupation amplitudes u k and v k have been shown in Fig. 
([2]). We note that in the limit of small interaction (V = 
and A = 0), the free energy (Eq. (fT2|) ) boils down to the 
energy of the noninteracting system of electron-hole pair 
states occupied up to the chemical potential (i.e., v k ^ and 

v k [ 2 equal to one below hw and zero above). 

Once the free energy is computed, the magnetization 
of the Eu moments is calculated within the mean-field 
approximation. Each Eu moment experiences an effective 
magnetic field H e / f given by 



E 



eff 



dF {M)/N 
DM 



(18) 
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FIG. 2: Electron-hole pair amplitudes u k ] and v k -\ as a func- 
tion of momentum k. Pairs of spin up electrons and spin 
down holes are filled up to the Fermi momentum kp-\ with- 
out any interactions. The step function is in the limit of zero 
interactions (Vo = 0, A — > 0), while both the Coulomb inter- 
action and the coupling to light smear out the Fermi surface. 
Here parameters are: A = 0.1 eV, m* = 1.5m e , Vo — 20 
eV.A 3 and Yiuj - E g = 0.33 eV, J M = 2.19 x 10" 5 eV and 
J = —0.036 eV . The pair amplitudes for the other spin, u k y 
and Vki, (not shown in the figure) are slightly shifted to the 
left, since fe_pj. is slightly lower than 



where the Bj is the Brillouin function and (3 — (kpT) -1 
and M(S) — g^sS. The above two equations are solved 
sclf-consistently for a given temperature using the free 
energy expression Eq. (fT2"]l . 

Note that in our formulation, we have neglected the 
thermal excitation of the electrons from the valence to the 
conduction bands, because the number of such carriers 
is much smaller than the number of the photogenerated 
carriers. The latter number, obtained by integrating v ka 
over all momentum, is typically ~ 1CP 4 / A 2 , while the 
number of intrinsic thermal carriers is many orders of 
magnitude smaller^ when the semiconductor bandgap 
is as large as ~ leV, as is the case for EuS. 




so that the magnetization is given by the expression 



FIG. 3: Dependence of the magnetization M as a function 
of temperature T with (dashed line) and without (solid line) 
the incident light. Parameters are: A = O.leV, m* = 1.5, 
Vo = 20 eV.A 3 , huj-E g = 1.5eV, J M = 2.19 x 10" 5 eV, and 
J = -0.036 eV. 



M(T) = g^ B SBj(pg^ B SH eff (M)), 



(19) 



The calculated magnetization and the number of pho- 
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FIG. 4: Pair occupancies as a function of temperature, n-f 
denotes the number of electron-hole pairs with spin up elec- 
trons and spin down holes, while nj. denotes the number of 
pairs with opposite spins. Below the critical temperature, n-f 
is larger than nj. because of the lower Zeeman energy of the 
former pairs. Above T c , there is no net megnetization, so that 
the energy of the electron is independent of the spin (zero Zee- 
man splitting), so that the pair occupancies are independent 
of the spin of the electron or hole. Parameters are the same 
as for Fig. 




FIG. 5: Photoinduced increase of the magnetic transition 
temerature as a function of the photon energy. Parameters 
are the same as in Fig. 2. A larger photon energy leads to a 
larger number of carriers which add to the ferromagnet 

toexcited electron-hole pair as a function of the temper- 
ature are shown in the Figs. [3J and 2] respectively. These 
are the two main results of this paper. As seen from 
the figures, there is a transition from the paramagnetic 
to the ferromagnetic state as temperature is decreased. 
Fig. [3] shows that with the incident light and for typ- 
ical parameters, the photoinduced increase in the tran- 
sition temperature T c is of the order of a fraction of a 
degree, AT C ~ 0.2K, in line with the experimental value 8 



AT C - 0.16K for the light frequency oihw-E g = Q.iieV. 

Fig. U shows the number of pairs of spin-up electron 
and spin-down hole and pairs of spin-down electron and 
spin-up holes as function of temperature. While the num- 
ber of the electron-hole pairs is controlled by the chemical 
potential (/i = fno) 1 the distribution of the pairs between 
the two different spin types (spin-up electron and spin- 
down hole or vice versa) is determined by the magnitude 
of (5 1 ), which affects the pair energy as indicated from 
Eqs. (fT"5f and (fl9|) . The effective magnetic field expe- 
rienced by the Eu spins is roughly proportional to the 
difference of the number of the two spin types, as may be 
seen by keeping the first line in the free energy Eq. (fT!?|) 
and using Eqs. (fTB"]) and (fT8]) , so that H e ff oc (n^ — nj). 
This in turn implies that this number difference and the 
magnetization M follow one another, both ultimately go- 
ing to zero at the transition temperature T c as seen from 
Figs. [3jand[4] In the paramagnetic phase, the photoex- 
cited electron- hole pairs are present, but the two spin 
types are the same in number. 

In Fig. [5] we plot the increase in T c vs. the inci- 
dent light frequency, which as discussed earlier effectively 
serves as the chemical potential. We see that as the fre- 
quency increases, more and more photocarriers are cre- 
ated and hence AT C increases. The figure shows that 
even when hu; — E g < 0, there is an increase in T c . This 
means that we have photoinduced magnetization even for 
the sub-bandgap light frequency, an issue that has been 
studied in detail by previous authors i 16 i 24 Note that we 
have neglected the relaxation of the excited carriers, so 
that in our case carriers are filled up to the chemical po- 
tential [i — hw at equilibrium. This relaxation effect will 
diminish the total number of photoexcited carriers, re- 
sulting in a somewhat lower value of AT C than predicted 
from our analysis. Such relaxation effects must be taken 
into account if a more quantitative description is desired. 

In summary, we studied the photoinduced enhance- 
ment of magnetization in the ferromagnetic semiconduc- 
tors by solving a model Hamiltonian which takes into 
account the Coulomb interaction, light-matter interac- 
tion, the exchange interaction between the carriers and 
the fixed magnetic moments, as well as the direct ex- 
change interaction among the localized moments. From 
the calculated free energy, we computed the magnetiza- 
tion of the system as a function of temperature. There is 
an enhancement of T c by a fraction of a degree K due to 
carriers created by light, in semi-quantitative agreement 
with the observed experimental results for EuS. 

This work was supported by the U. S. Department of 
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